Gas condensation within a bundle of carbon 
nanotubes — effects of screening. 



M.K.Kostov, J.C.Lewis and M.W.Cole 

Department of Physics, 
The Pennsylvania State University, University Park, PA 16802, USA 

February 1, 2008 



One of the principal themes of condensed matter physics is the sensitivity of a system's behavior 
to its dimensionality |Q, H |^j. Much interest in surface physics is attributable to the rich diversity of 
possible two dimensional (2D) phase diagrams. The discovery of the carbon nanotube (NT) has created 
the possibilty of interesting ID and quasi - ID phases involving adsorbed atoms and molecules [4-9]. In 
the case of small molecules/atoms, such as H2, He and Ne, the most energetically favored site is believed 
to be the interstitial channel (IC) between NTs (To), which is indeed a ID domain. In this paper we 
evaluate the phase behavior of a system of many H2 molecules within the ICs. An important role is 
played by the carbon environment, which screens the attraction between hydrogens. As a result, we find 
that the ground state cohesive energy of the system is reduced to a very small value ( less than 5% of 
the well depth of the free space pair potential ). In this paper, we estimate the critical temperature 
for the condensation of hydrogen using alternatively localized and delocalized models of the hydrogen. 
A comparison is presented with similar results for the adsorption of neon and helium. In all cases, the 
effects of screening in reducing the transition temperature are significant. 

Our study is a successor to a recent investigation of the intrachannel screening j?J . That project eval- 
uated the Axilrod-Teller-Muto triple dipole (DDD) interaction between two hydrogen molecules within 
the same IC and the surrounding carbon atoms. The DDD interaction has been shown empirically to 
account quantitatively for the three-body interaction involving simple gases in 3D and semiquantitativcly 
with similar effects in adsorbed films Jll], The resulting net interaction is: 

V total (r)=V(r) + V DDD {r), (1) 

where V(r) is the free space interaction. Vtotai{r) has a well depth, 1.35 meV, which is just 46% of the 
free space well depth. As a consequence of the reduced attraction the H2 dimer's binding energy falls 
from 4.5 K to zero. This absence of the dimer bound state implies that the ground state of the strictly 
ID system is a gas || . This result is manifested in the absence of a finite p minimum in the ID equation 
of state Em(p) at temperature T — in Fig. 1 (dashed curve); here p is the ID density. This curve 
has been obtained by Boronat and Gordillo jl3| using their diffusion Monte Carlo method, together with 
the screened interaction of Kostov et at (t]] . The curve contrasts strikingly with the result of a previous 
study which found a liquid with cohesive energy 4.83 K for the ID system and the free space interaction 
between H 2 molecules []. 

To study the collective properties of the H 2 , we must take into account the interaction between 
molecules of neighboring channels. Fig. 2 depicts both the unscreened and the screened intermolecular 
interactions, the latter including the DDD term which we evaluated using the same procedure as Kostov 
et al. JtJ. Similarly, Fig. 3 and Fig. 4 illustrate the DDD screening effect in the interchannel interaction 
for the case of He and Ne atoms, adsorbed within the ICs. In all three cases, we observe quite a large 
DDD effect, which alters significantly the free space interchannel potential. Note that the distance scale 

1 Because of some uncertainty about the adequacy of the DDD in describing the screening, we have explored the conse- 
quence of artificially reducing this interaction. Multiplying this interaction by 0.845 ± 0.01 yields the threshold for there to 
be a dimer bound state. Multiplication by one-half yields a dimer binding energy of 1.1 K . 



in Fig. 2, Fig. 3 and Fig. 4 is the longitudinal (along the axis of the IC) displacement z between the two 
particles of neighboring channels. The spacing between neighboring ICs is a fixed parameter d — 9.8 A. 

One can easily obtain a rigorous upper bound for the interacting system's energy. To do so, we identify 
the system of uncoupled H2 (i.e. independent ICs) as the unperturbed system. The ground state of this 
latter system is a set of uncondensed ID gases. The variational theorem asserts that the fully interacting 
system's energy satisfies 

E{p) < E mr {p) , (2) 

E var (p) = E lD {p) + ^ [°° dzV total \{z 2 + p 2 ) 1 ' 2 ] . (3) 
* J-00 L 

The prefactor of 3/2 is derived from a coordination number of 3 neighboring ICs and the need to 
avoid double counting. More distant tubes (omitted here) would further lower this variational energy by 
less than 1%, which we neglect. Note that in Fig. 2 the DDD term has regions which are both positive 
and negative. The integral is such that the negative contribution outweighs the positive contribution. 
This means that the environment enhances the interaction between molecules in adjacent channels. This 
might be called an "antiscreening effect" . 

Fig. 1 presents the resulting variational energy. There is a minimum at density po ~ 0.091 A -1 , at 
which the ground state energy per particle is 0.15 ±0.01 K . This is a remarkably low density and weakly 
bound liquid. It is analogous to that found recently for 4 He in the same environment The latter 

calculation does not include environmental screening however. 

What is the phase behavior of this unusual system? As was found for 4 He, there is a low temperature 
region of separation into liquid and vapor phases. A droplet of "liquid" will involve condensation in 
contiguous ICs; this phase will be strongly anisotropic, and for that reason alone is worth further study. 
The critical region should be 3D Ising-like, since the transition is inherently 3D. The critical temperature 
may be estimated in at least two ways, based on alternative assumptions. The first of these ways is 
simply to note the general similarity between cohesive energy and critical temperature. For 4 He in 3D 
for example, the cohesive energy is E c = 7.17 K and the critical temperature is 5.2 K. In 2D, these values 
are also similar (E c = 0.85 K, T c = 0.87 K)jy|. On this crude basis, we estimate T c ~ 0.1 K for H2 in 
NTs. 

However, a quite different estimate is obtained if the H2 molecules become localized at sites created by 
the periodic potential provided by the NTs. In this case, a lattice gas model is appropriate. To describe 
the phase behavior of H2 in this model we follow the procedure of Cole et al. H]. In that approach 
we obtain J t w 14 mK and c w 0.007, where an anistropic simple cubic Ising model was assumed with 
an interaction strength J z between neighboring spins along the z axis (within the same channel) and 
a transverse interaction J t = cJ z . The critical temperature of the H2 condensation transition in this 
lattice-gas model is then T c ~ 1.19 K p6| . 

Analogous assessments of the condensation properties of He and Ne in that model yield Jt ~ 2.6 mK, 
c « 0.009 and T c ~ 0.18 K for He , and J t ~ 15 mK, c « 0.009 and T c ~ 1.05 K for Ne. The previous 
study of He condensation found T c ~ 0.36 K in the nonscreening case. 

In summary, we have derived estimates of the transition temperatures for quasi- ID condensations 
of gases in nanotube bundles. The effect of the carbon environment is significant. It provides strong 
screening of the mirachannel attraction but overall enhances the mterehannel attraction n. These phase 
transitions would be dramatically observable in specific heat measurements at low temperature because 
the tubes are thermodynamically inert in that temperature regime. 
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2 If the screening of both intra- and interchannel interactions were omitted, we would find the following values of T c for 
H2: ~ 4.8 K in the quasifree model and ~ 3.01 K in the localized model. 
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Figure Captions 



Fig. 1. Energy per molecule of ID H 2 within a single channel (dashed curve from Ref.[jl3|) compared 
with the energy per molecule in the case when interchannel interaction is included (full curve). 

Fig. 2. Potential energy of interaction between two H2 molecules in free space (dashed curve) and 
in adjacent channels (full curve) are shown. The difference arises from the DDD interaction of the two 
H2 molecules and the surrounding carbon atoms (dotted curve). The abscissa is the difference in z 
coordinates of the two molecules. 

Fig. 3. Same as Fig. 2 but for He. 

Fig. 4. Same as Fig. 2 but for Ne. 
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